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ABSTRACT 

Simulations of an isolated Milky Way-like galaxy, in which supernovae power a galactic 
fountain, reproduce the observed velocity and 21 cm brightness statistics of galactic 
neutral hydrogen (Hi). The simulated galaxy consists of a thin Hi disk, similar in 
extent and brightness to that observed in the Milky Way, and extra-planar neutral gas 
at a range of velocities due to the galactic fountain. Mock observations of the neutral 
gas resemble the Hi flux measurements from the Leiden- Argentine-Bonn (LAB) Hi 
survey, including a high-velocity tail which matches well with observations of high- 
velocity clouds. The simulated high- velocity clouds are typically found close to the 
galactic disk, with a typical line-of-sight distance of 13 kpc from observers on the 
solar circle. The fountain efficiently cycles matter from the centre of the galaxy to its 
outskirts at a rate of around O.5M0yr~^. 
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1 INTRODUCTION 

Radio observations of the 21 cm hydrogen emission line re- 
veal that the Milky Way (MW) contains a thin Hi disk 
surrounded by a population of 'clouds' with velocities in- 
compatible with models of galactic rotation yet apparently 
. not part of the Hubble flow either (Muller et al (1963)). The 
nature of these 'high- velocity clouds' (HVCs) remains some- 
, what unclear mostly because it is difficult to determine their 
■ distances and hence infer physical properties. 

Observations of stars with known distances along the 
line of sight to a HVC can be used to constrain the distance 
to the cloud, by testing whether or not it is detected in 
absorp tion in the stellar spectrum (see e.g. ISchwarz et all 
(|1995|)). Unfortunately such constraints are avai l able f or 
only a relative ly small number of HVCs (Wakker| (|200lh ). 
IPutman et aD (|2003) use constraints from Ha emission, and 
find that most HVCs are within 40 k pc, except fo r those 
associated with the Magellanic stream. TSruns et"al] (|200lh 
claim that some fraction of HVCs display a head-tail mor- 
phology, may be a result of inte r action with a diffuse ambient 
galactic wind (|Quilis k, Moord (|200lh l. Searches for stars as- 
socia ted with HVCs ha ve so far resulted in non-detections 
(e.g. IHodp et all (|200/t )). Metallicities have been measured 
for a small numbe r of H VCs and vary over a wide range 
(e.g. iGibson et all (|200lh l. suggesting that HVCs are not 
a homogeneous set. Finally one can observe other galax- 
ies to infer the nature of HVCs from their projected dis- 
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tances. The Andromeda galaxy (M31) ha s a popu l ation of 
HVCs close (^ 50 kpc) to its disk (Westmeier et al.1 (|200/t )). 
Other nearby sp iral galaxies al s o con ta in extra-planar neu - 
tral gas (see e.g. iBarbieri et al. I (|2005l ll. |Pisano et al.1 (|2007l ) 
have searched for HVC analogues in six loose groups of 
galaxies, similar to the Local Group. Their failure to de- 
tect compact HVCs implies that any H i clouds are near 
to the galaxies in this group: they are not roaming freely 
throughout the group itself. These observations do not un- 
ambiguously determine the nature of HVCs, and several the- 
oretical models for them have been proposed. Oort (1970) 
discusses a model where the HVCs are close to the MW disk, 
and r esult from of a 'galactic fountain' (|Shapiro Fieldl 
(| 19761 ) V Gas from the galactic disk rises buoyantly after 
being heated by supernovae (SNe), becomes thermally un- 
stable and cools radiatively into neutral clouds. Once the 
clouds are dense they fall back ballistically onto the disk, 
and are seen as the HVCs. Whether SNe explosions can 
start a galactic fountain has been investigated both theo- 
retically (.Kahn, (^1998.)) and through numeri cal simulations 
(|de Avilled (|l998h : lde Avillez k Berrvl (|200lh ). 

Blitz et al. (1999; see also lBraun k BurtonI (|l999l )) sug- 
gest (some) HVCs are neutral gas associated with the nu- 
merous dark matter substructures seen in simulations of 
haloes of galaxies and groups of galaxies. Such HVCs are 
at large distances (> 40 kpc) from the g alactic centre, 
and d istributed throughout the local group. IConnors et all 
(|2006l l report the detection of Hi clouds in fully cosmo- 
logical ACDM simulations of galaxy for mation, and their 
HVCs are reminiscent of Blitz's clouds. ^ Mailer k Bullockl 
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discuss how a population of clouds at distances of 
~ 150 kpc may result from the halo gas becoming ther- 
mally unstable. Given the rich variety HVC properties it 
seems likely that more than one of the above models, with an 
additional component resulting from tidal interactions (e.g. 
iFerguson et all (|2006l ll is required to explain the full cloud 
distribution. In this letter we perform numerical simulations 
of star formation in an isolated Milky Way-like galaxy, in 
which SNe power a galactic fountain. We analyse the sim- 
ulation in terms of mock Hi observations and demonstrate 
they reproduce many of the observed features of the H i disk 
and its HVCs. 

This letter is structured as follows. In section [2] we pro- 
vide details of the numerical scheme used in our investiga- 
tion and in section [3] we present a comparison between our 
simulations and the distribution of galactic Hi. 



2 NUMERICAL METHOD 

2.1 Star Formation and Feedback model 

The star formation and feedback prescription s used in the 
simulations is described in detail inl Booth et all (|200/t ) (here- 
after BTO07), here we briefly review its main features. The 
scheme treats the interstellar medium (ISM) in terms of 
three distinct but interacting components: cold (T < 10^ K) 
molecular clouds surrounded by a warm (T ^ 10^ K) am- 
bient phase in the disk, interspersed by a hot (T ~ 10^ K) 
tenuous phase that extends into the halo, and is powered 
by SNe. Gas is cycled through these phases by a number of 
processes. 

Conversion of gas to stars is regulated by the rate at 
which molecular clouds form and get destroyed. Clouds 
form from thermally unstable ambient gas and get de- 
stroyed directly by feedback from massive stars, and indi- 
rectly through thermal conduction. Gas cooling is due to 
Compton and line coolin g, using inte r polat ion tables gen- 
erated using Cloud^Q ( Per la nd et all (|l998h V Feedback in 
the ambient gas phase cycles gas into a hot galactic fountain 
or wind. We model the ambient and hot gas phases hydro- 
dynamically using smoothed particle hydr odynamics (SPH, 
iGingold MonaghanI (|l977l ): ILucvI (1 19771)1 as impl emented 
in the Gadget2 simulation code (Springel' ("2005')). Moti- 
vated by the fact that we cannot resolve the Jeans mass 
of the molecular gas, the clouds are modelled using 'sticky 
particles' that move ballistically through the ambient gas, 
but may coagulate when colliding. When such a cloud has 
built-up enough mass through accretion and coagulation to 
form a Giant Molecular Cloud (CMC), it collapses into stars, 
which then destroy the CMC through stellar winds and SNe 
explosions. 

BTO07 show that this model produces a multiphase 
medium with cold clouds, a warm disk, hot SN bubbles and 
a hot, tenuous halo, similar to that observed in spiral galax- 
ies. The star formation rate, surface Hi density, molecular 
fraction, and molecular cloud mass-spectrum of simulated 
galaxies match closely those observed in the MW. 



^ We thank our colleagues J Schaye, C Dalla Vecchia and R 
Wiersma for allowing us to use these routines. 



2.2 Model Galaxy 

The initial conditions for our MW-like galaxy are gener- 
ated using a publicly ava ilable programme, GalactICs 
('Kuiiken Dubinskil (|l995l )). which generates a near equi- 
librium galaxy consisting of an approximately exponential 
disk, a bulge and a (dark matter) halo, using collision- 
less particles. The bulge, disk and halo have masses in 
the ratio 0.31:1.00:28.27 (giving a baryon fraction of ^0.3 
times the universal value). The total mass of the system is 
2.21 X lO^^M© and the circular velocity at the solar radius, 
r = 8.5 kpc, is - 220km s"\ 

We convert 10% of the disk particles into SPH parti- 
cles with temperature lO^iC, the remaining 90% of the disk 
represents a stellar disk with population ages set by assum- 
ing that the galaxy had a constant star formation rate of 
lM0yr~^. The bulge is changed in its entirety into an old 
stellar population. Finally 1% of the material in the halo is 
changed into SPH particles with a temperature of l(fK. The 
galaxy is now no longer in equilibrium and starts converting 
atomic gas into molecular clouds, which in turn form stars, 
yet soon settles into a near-equilibrium state where feedback 
regulates the star formation rate at approximately lM0/yr. 
It was shown in BTO07 that the resulting galaxy resem- 
bles the MW in terms of its gas fraction, molecular fraction, 
gas distribution and Schmidt-Kennicutt star formation law. 
Our highest resolution simulation initially contains ^ 3x 10^ 
particles in the disk (particle mass 1.5 x lO^M©), but it 
has been demonstrated in BTO07 that that the properties 
of this galaxy are not strongly dependent on resolution. 



2.3 Mock Hi observations 

We assume that the Hi fraction of the ambient gas is set 
by the balance between photo-ionisation, collisional ionisa- 
tion, and recombinat ions, which we compute using Cloudy 
(|Ferland et al.l (|l998l)), imposing the UV back ground at red- 
shift z = given bv lHaardt &: Madaul (l200ll ) . This Hi gas 
can be detected by its hyperfine emission line ( Van de HulstI 
(1945)). A hydrogen atom in its ground state has higher 
energy when the spins of proton and electron, Jp and Je, 
are parallel than when they are anti-parallel. Spontaneous 
transitions to the lower energy state occur at a rate gov- 
erned by the Einstein coefficient, Aiq = 2.85 x 10~^^ s~^, 
and result in the emission of a '21 cm photon' of wave- 
length A = 21.11cm = c/v. For the densities, nn, and tem- 
peratures, T, relevant to the ISM, the typical time r be- 
tween collisions, r - 100 x (riH/cm"^)"^ (T/IO'^K)"^/ V, 
is much shorter than the spontaneous 21cm transition rate, 
Til 2 — s ~ 10^ yr, and so collisions keep the levels of the 
low and high states, no and ni, in the ratio ni = 3no. The 
mean emission rate of 21 cm photons per neutral hydrogen 
is then given by Ai^m = (3/4) Aio. The 21 cm flux of a 
cloud with neutral hydrogen mass Mhi at distance D is 

(3/4) Aio (Mm/mH) 

see e.g. ISpitzerl (|l978h ). In radio astronomy, this is quoted 
in terms of the line flux, f S{iy) dv, with S{iy) expressed 
in Jy and the line- width in km s~^, F — J S{iy)diy = 
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f S{iy) (y jc) dv, hence (|Wakker van WoerdenI (|l99lh ^ 



y L Jy J Lkms-iJ 



Jy kms ^ V 

1 Mhi /kpc\2 
0.235 Mq\D) ■ 



(2) 
(3) 



The 21cm 'brightness temperature' ,Ts, of an object is de- 
fined as the temperature at which a black-body emits the 
same flux. The conversion from flux to brightness temper- 
ature is then given by Tb/S — R, where the telescope- 
dependent conversion factor we use is = 0.158K Jy~^ in 
order to match the observational survery we are comparing 
our results against (|Hulsbosch Wakkerl (p^988)). To com- 
pute the simulated H I flux we place an observer in the sim- 
ulated galaxy and evaluate the net flux received by an ideal 
radio telescope with a beam size 0. For a single SPH particle 
at position r^, the fraction dm/m of mass that falls within 
the beam at distance between r and r + is 

dm/m^ (^27rr^(cos(6')-l)(ir^M^(|r^-r|,/i), (4) 

where h is the smoothing length of the particle, and W 
the SPH kernel. The total flux received from this particle 
is computed from Eq. ([T|), and is represented by a Gaus- 
sian emission line centred at velocity v • r/r with width 
(Jv — (k^T /rrihY^'^ . The total spectrum is obtained inte- 
grating over dr and summing over all particles. 

Simulations have been performed at two different mass 
resolutions to ensure that resolution effects do not affect 
our results. Observations have been repeated for a number 
of observers along the solar circle, and at times spanning 
a period of IGyr. We use these observers to compute er- 
ror bars on the mock observations. The distribution of H i 
at galactic latitudes, 6, greater than 20° does not depend 
strongly on the time at which observations are made, show- 
ing that the galaxy our results do not represent a transient 
feature. Within the galactic disk (h < 20°) the mean bright- 
ness temperature of the gas decreases slowly over time as 
the gas disk is converted into stars. 



3 RESULTS AND DISCUSSION 

The all-sky 21 cm brightness distribution of the simulated 
galaxy looks remarkably similar to the observed H i in the 
Milky Way, as measured by the Leiden- Argentine-Bonn 
fLAB. lKalberla at all (|2005l )) Hi survey (Fig.[l]). The LAB 
survey has angular resolution of 0.5° x0.5° but unfortunately 
our numerical simulation does not have enough particles to 
resolve structures on such small scales. We calculate the 
mean angular extent of the particles that contribute to the 
flux in a given direction on the sky in the simulated galaxy, 
and then smooth the LAB survey with a Gaussian kernel 
of the same width. At low galactic latitudes this smoothing 
angle is typically less than 1°. However, at high latitudes 
the mean smoothing length is much larger, and at \b\ > 60° 
can reach up to 20° due to the relatively small number of 
SPH particles at these latitudes. 

Both observed and simulated brightness maps display a 
bright and thin H I disk in the plane of the MW, embedded in 
a thicker cooler envelope (Tb ~ 10^ K), with an even cooler 
component (Tb ^ 100 K) at high velocities, It'isrl > 100 km 
s~^, with respect to the local standard of rest (LSR). The 
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Figure 2. Velocities with respect to the local standard of rest of 
Hi gas in the LAB survey (top panel) and the simulated galaxy 
(bottom panel) as function of galactic longitude, /. Intensity on 
this plot is the 21cm brightness of the H I , integrated over galac- 
tic latitude. The simulated galaxy has H i gas out to velocities of 
\v\sj-\ ~ 200 km s"-*^ in a characteristic pattern also seen in the 
LAB. The blue points (bottom panel) are HVCs from the cata- 
logue of Lockman et al (2002). The simulation does not resolve 
small clouds, but the velocity distribution of the simulated Hi 
traces the same regions in (/, vi^j.) space as the observed clouds. 

brightness temperature Tb, and its fall-off with latitude, is 
very similar in the observed and simulated maps. The min- 
imum H I brightness temperature in the simulated map is 
95K, in good agreement with Hi observations of the MW, 
where it is found that every line-of-sight contains easily ob- 
servable H I. The simulated high velocity gas (It'isrl > 100 km 
s~^) forms a nearly uniform background. Due to the rela- 
tively small number of particles that are flagged as high 
velocity at any one time the spatial resolution is poor, espe- 
cially at high latitudes, and these simulations do not resolve 
the flne structure seen in the LAB. When the LAB data are 
smoothed to the same resolution as the simulation the re- 
sulting distribution matches closely, with a mean brightness 
temperature of Tb ~ 23.22 K, as compared to 20.33 K in 
the simulations. 

The velocity distribution of the simulated Hi also 
matches well with the LAB data (Fig [2]). Although we can- 
not resolve individual HVCs, the simulated Hi velocities 
are in go od agreement with th e properties of the HV C cat- 
alogue of ' Lockman et all (|2002l l . [Lockman et~al ^00^ iden- 
tifled some of the detections in this survey with external 
galaxies, we have removed these from the plot. Additionally, 
clouds that were identified as being part of the Magellanic 
stream, which dominate the extreme negative velocity flow 
(|Mathewson, Clearv, Murravl (| 19741 ) 1, were removed. 

The line flux, J S{iy)diy, for high velocity gas, \v\sr\ > 
lOOkms"^, is shown in Fig. O Solid black lines represent 
the results of I Wakked i 199lh with the emission due to the 
Magellanic stream and outer arm of the MW removed, as 
we do not expect an isolated galaxy to match these fea- 
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Velocity Integrated HVCs 




Figure 1. All-sky maps of 21 cm brightness temperatures in the LAB survey (top panels) and in simulations (lower panels), integrated 
over all velocities (left panels) or integrated over high velocities {\v\q^\ > 100 km s"-*^) only (right panels). The observational data has 
been convolved with a spatially adaptive Gaussian, to make the angular resolution of both the simulated and observed data the same. 
Both observed and simulated galaxy consist of a thin (6 < 5°) Hi disk with brightness temperature Tb > 5000 K, embedded in a thicker 
(6 ~ 45°) disk with Tb ~ 1000 K, with Tb ~ 100 K at higher latitudes. The high velocity gas has Tb < 100 K in both observed and 
simulated galaxy. The simulation does not reproduce the smaller-scale structure seen in the LAB due to lack of numerical resolution. 




Figure 3. Line flux of 21 cm emission in galactic longitude I (left panel) and latitude b (right panel) for high velocity gas. The solid black 
lines are from the survey for HVCs by Wakker (1991). The observational data has the populations of clouds belonging to the Magellanic 
stream and outer arm removed, as an isolated MW galaxy can not reproduce these. The grey histogram is obtained from the simulations, 
with error bars a measure of the scatter between simulated observers at the solar circle. The simulated galaxy looks very similar to the 
observed one in I, but the distribution in b is slightly more extended in the real MW, although the simulated scatter is large. 



tures. Errors on the mock data are calculated by repeating 
the measurements for observers at different points along the 
solar circle. The mock and real data look remarkably similar. 
The distribution of neutral gas perpendicular to the galac- 
tic plane is approximately exponential with a scale he i ght o f 
5 kpc, in agreement with the predictions o f iBregmanI (jl980l l 
for a galactic fountain. Although as noted in Fig [3] the distri- 
bution of 21cm emission in galactic latitude is slightly more 
concentrated in the simulated data than the observed data, 
when averaged over the whole sky the mean values of the 
smoothed maps agree to within 10% and are 485K for the 
observed data, and 446K for the simulated data. 

It is difficult to measure the distances to HVCs in the 
real universe. However, in our simulations this information 



is preserved and can be measured easily (Fig|4|. 50% of the 
flux is emitted within a distance of 13 kpc from our ob- 
serve r. This figure i s in li n e with that observed other galax- 
ies ( Barbieri etHI (l2005l ) : IPisano et a l. (2007)), notably in 
M 31 bv lWestmeier et all (|2005l l. who found that HVCs were 



generally at a projected distance of less than 15 kpc. 

Gas in the fountain rains back on the disk at a radius 
(^returned) which is generally larger than the radius (rejected) 
from where it was launched (Fig[5|. An SPH particle is clas- 
sified as being part of the fountain if it is ejected to a vertical 
distance of more than 2 kpc from the galactic plane and at 
a later time falls back within 1 kpc of the galactic disk. The 
galactic fountain has an overall effect of moving gas from the 
inner part of the galactic disk to its outskirts: most fountain 
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Figure 4. Flux weighted mean distance to H I emitting high 
velocity (l-uisrl > 100 km s""*^) gas in the simulated galaxy. Fifty 
per cent of the emission occurs within ~ 13 kpc (vertical line). 
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Figure 5. Histogram showing the relation between the galac- 
tic radius at which gas is ejected from the galaxy (rejected) 5 
and the radius at which it next passes through the galactic disk 
(^returned)- Overall the galactic fountain acts to move gas from 
the centre of the galactic disk to its outskirts. 

particles with rejected < 8 kpc have rr etumed/^eiected ^, 4 

hence rain back outside the solar circle. ICorbelli Salpeterl 
(Il988l) note that the effect of this 'outer galactic fountain' 
depends upon the physical properties of the gas that returns 
to the disk at large radii and could either evaporate the H i 
disk, or cause it to grow. The rate at which gas returns to the 
outer disk in the simulations, M ^ 0.5Mq yr~^, is in agree- 
ment with observational estimates and could b e enough to 
drive the observed turbulence (Santilla n et al.l f2007)). We 
find that a peak mass flux of at least 0.5Mq yr~^ through 
the disk out to a radius of 15kpc, suggesting that a galactic 
fountain of this form is capable of driving ISM turbulence all 
the way to the outskirts of the galactic disk. In our simula- 
tions, M decreases gradually as the quiescent star formation 
rate gradually decreases due to the depletion of the local gas 
supply. However, we note that in a fully cosmological set- 
ting bursts of star formation due to galaxy mergers coupled 
with infall of gas from the surrounding intergalactic medium 
could continue to drive this process over the lifetime of the 
galaxy. 

We conclude that our multi-phase star formation imple- 
mentation naturally produces a galactic fountain in a Milky 
Way-like model galaxy. The neutral hydrogen in the foun- 



tain has a spatial and velocity distribution in good agree- 
ment with a variety of observations for the MW and other 
spiral galaxies. 
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